QED 3 on a space-time lattice: compact versus noncompact formulation 
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We study quantum electrodynamics in a (2+l)-dimensional space-time with two flavors of dy- 
namical fermions by numerical simulations on the lattice. We discretize the theory using both the 
compact and the noncompact formulations and analyze the behavior of the chiral condensate and of 
the monopole density in the finite lattice regime as well as in the continuum limit. By comparing the 
results obtained with the two approaches, we draw some conclusions about the possible equivalence 
of the two lattice formulations in the continuum limit. 
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I. INTRODUCTION 

Quantum electrodynamics in 2+1 dimensions (QED3) 
is interesting as a toy model for investigating the mech- 
anism of confinement in gauge theories p|, and as an 
effective description of low-dimensional, correlated, elec- 
tronic condensed matter systems, like spin systems - 
or high-T c superconductors Q. While the compact for- 
mulation of QED 3 appears to be more suitable for study- 
ing the mechanism of confinement, both compact p| and 
noncompact formulations arise in condensed matter sys- 
tems. Our paper aims to elucidate some aspects of the 
relationship between these two formulations of QED3 on 
the lattice. 

Polyakov showed that compact QED 3 without fermion 
degrees of freedom is always confining 0. Any pair of 
test electric charge and anti-charge is confined by a lin- 
ear potential, as an effect of proliferation of instantons, 
which are magnetic monopole solutions in three dimen- 
sions. The plasma of such monopoles is what is respon- 
sible for confinement of electrically charged particles. If 
compact QED3 is coupled to matter fields it has been 
argued @ that the interaction between monopoles could 
turn from 1/x to —In (a;) at large distances x, so that 
the deconfined phase may become stable at low tem- 
perature. The issue of the existence of a confinement- 
deconfinement transition in QED3 at T = is still con- 
troversial, as it has also been proposed that compact 
QED3 with massless fermions is always in the confined 
phase 0, @ ; also, in the limit of large flavor number, it 
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has been argued that monopoles should not play any role 
in the confinement mechanism 9] . At finite temperature, 
parity invariant QED 3 coupled with fermionic matter un- 
dergoes a Bcrczinsky-Kosterlitz-Thouless transition to a 
deconfined phase [lCj . 

The issue of charge confinement in 2 + 1 dimensional 
gauge models comes out to be relevant in the context 
of quantum phase transitions, as well. Indeed, recently 
it has been proposed that phenomena similar to decon- 
finement in high energy physics might appear in planar 
correlated systems, driven to a quantum (that is, zero- 
temperature) phase transition between an antiferromag- 
netically ordered (Need) phase, and a phase with no order 
by continuous symmetry breaking 0, • The most suit- 
able candidate for a theoretical description of the system 
near the quantum critical point is a planar gauge the- 
ory, either with Fermionic matter 0, or with Bosonic 
matter 

At finite T noncompact QED3 comes about to be rel- 
evant in the analysis of the pseudo gap phase °f 
cuprates. This phase arises from the fact that, upon 
doping the cuprate, a gap opens at some temperature T* 
which is quite larger than the critical temperature Tc 
for the onset of superconductivity. Both temperatures 
T* and Tc are doping dependent quantities and the gap 
is strongly dependent upon the direction in momentum 
space, since it exhibits <i-wave symmetry [12J. 

In Fig. ^ we report the phase diagram of high-T c 
cuprates. For small- a; phase is characterized by an 
insulating antiferromagnet (AF); by increasing x, this 
phase evolves into a spin density wave (SDW), that is 
a weak antiferromagnet. The pseudo gap phase is lo- 
cated between this phase and the d-wave superconduct- 
ing (dSC) one. 

The effective theory of the pseudo gap phase turns 
out to be QED 3 0, [lj, [l5| , with spatial anisotropies in 
the covariant derivatives, that is with different values for 
the Fermi and the Gap velocities , and with Fermionic 



matter given by spin-1/2 chargeless excitations of the 
superconducting state (spinons). These excitations are 
minimally coupled to a massless gauge field, which arises 
from the fluctuating topological defects in the supercon- 
ducting phase. The SDW order parameter is identified 
with the order parameter for chiral symmetry breaking 
(CSB) in the gauge theory, that is, (^ji/j) E3- There can 
be two possibilities; if (iptp) is different from zero, then 
the d-wave superconducting phase is connected to the 
spin density wave one (see Fig. Q case b) ; otherwise the 
two phases are separated at T = by the pseudo gap 
phase (see Fig.^case a). 

Confinement and chiral symmetry breaking go essen- 
tially together as strong coupling phenomena in gauge 
theories; while confinement is an observed property of 
the strong interactions and it is an unproven, but widely 
believed feature of non-abelian gauge theories in four 
space-time dimensions, chiral symmetry is only an ap- 
proximate symmetry of particle physics, since the up and 
down quarks are light but not massless. Central to our 
understanding of CSB is the existence of a critical cou- 
pling: when fermions have a sufficiently strong attrac- 
tive interaction there is a pairing instability and the en- 
suing condensate breaks some of the flavor symmetries, 
generate quark masses, and represents chiral symmetry 
in the Nambu-Goldstone mode [Til Il7j . The issue of 
a critical coupling has been widely investigated in 2+1 
dimensional gauge theories [II [HHi. Typically, the 
dimensionless expansion parameter is 1/Nf. Using the 
Schwin ger- Dyson equations p"s| or a current algebra ap- 
proach [21j for QED3 and QCD3 one finds that there is 
a critical number of flavors, Nf C) such that only for Nf 
lesser than Nf tC chiral symmetry is broken; for Nt bigger 
than Nt c chirality is unbroken and quarks remain mass- 
less. For QEDg this result has been the subject of some 
debate [Il|llMll3j2lJMJli; there are, however, nu- 
merical simulations [271 l28l l29j of QED3 , which find an 
Nf iC remarkably close to the results reported in Ref. 0] . 

Even if far from the scaling regime, strong coupling 
gauge theories on the lattice provide interesting clues on 
the issue of CSB. In fact, one can show that, in the strong 
coupling limit, a Hamiltonian with N c colors of fermions 
and Nf/2 lattice flavors of staggered fermions is effec- 
tively a U(Nf/2) quantum antiferromagnet with repre- 
sentations determined by N c and Nf [3(|. CSB is then 
associated [3(| either to the formation of a U(l) commen- 
surate charge density wave or of a SU(Nf/2) spin den- 
sity wave, i.e. to the formation of Neel order. Quantum 
antifcrromagnets with the representations considered in 
Ref. have been analyzed in Ref. [3l| where it was 
found that, for small enough Nf, the ground state is or- 
dered. Also, when Nf is increased there is a phase transi- 
tion, for Nf ~ N c , to a disordered state. In this picture, 
the large N c limit is the classical limit where Neel order 
is favored and the small N c and large Nf limit are where 
fluctuations are large and disordered ground states are 
favored. 

We shall not try to ascertain in this paper the critical 
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FIG. 1: Phase diagram in the (x,T) plane [3; here x repre- 
sents the doping and T the temperature. 

number of flavours Nf tC . Here, we shall analyze the re- 
lationship between monopole density and fermion mass 
and compare the results obtained for the compact and 
noncompact lattice formulation of this gauge model. In 
particular, we revisit the analysis of Fiebig and Woloshyn 
of Refs. [33, H3| , where the dynamic equivalence between 
the two formulations of (isotropic) QED 3 is claimed to 
be valid in the finite lattice regime. In this paper we shall 
extend the comparison to the continuum limit, following 
the same approach as in Refs. j33, H3 : namely we shall 
analyze the behavior of the chiral condensate and of the 
monopole density as the continuum limit is reached. 

In Section II we describe the model and its properties 
both in the continuum and on the lattice. Moreover, the 
method for detecting monopoles on the lattice is illus- 
trated. 

In Section III a description of both compact and non- 
compact formulations of QED3 is given. 

In Section IV we present our numerical result for the 
chiral condensate and the monopole density in the region 
in which the continuum limit is reached. Then, we com- 
pare our results with those of Fiebig and Woloshyn [3^, 

Section V is devoted to conclusions. 



II. THE MODEL AND ITS PROPERTIES 

The continuum Lagrangian density describing QED3 
is given in Minkowski metric |34| by 

C = + AiD^i ~ mo&V< , (1) 

where Z)„ = dn — ieA^, F„ v is the field strength and 
the fermions ipi (i = 1, . . . , Nf) are 4-component spinors. 
Since QED3 is a super-renormalizable theory, dim[e] = 
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+1/2, the coupling does not display any energy depen- 
dence. One may define three 4x4 Dirac matrices 
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and two more matrices anticommuting with them: 
namely 



3 ,[01 
T = M i o 
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(3) 



The massless theory will therefore be invariant under the 
chiral transformations 



If one writes a 4-component spinor as ^ = 
mass term becomes 



?/>2 



(4) 
the 



Since in three dimensions the parity transformation reads 

lpl(x ,X\,X2) -> (Tilp2{XQ,-X\,X2) , 

ip2(xo,x 1 ,x 2 ) -> cti^i(x , -xi,x 2 ) , (5) 

then mipip is parity conserving. 

The lattice Euclidean action [2^, |35| using staggered 
fermion fields x> X; is given by 



JV 



S = S G + ^2^2xi(n)M n , m Xi(m) 



(6) 



j— 1 n,m 



where is the gauge field action and the fermion matrix 
is given by 

M n , m [U] (7) 
= E„=i, 2 ,3 2^ {[Uu(n)]6 m , n+ - [U£(m)]5 m , n -i>} ■ 

The action © allows to simulate N = 1 , 2 flavours of 
staggered fermions corresponding to = 2, 4 flavours 
of 4-component fermions tp j3||. Sq is different for the 
compact and noncompact formulation of QED3. 

For the compact formulation one has 



S G \U]=0 E 



1- 2 (^N + ^.N) 



(8) 



where U^ u (n) is the "plaquette variable" and — 
l/(e 2 a), a being the lattice spacing. Instead, in the non- 
compact formulation one has 



(9) 



where 

F^(n) = {a v (n+jl)-a v (n)}-{a^(n+v)~a tl (n)} (10) 

and a /i (n) is the phase of the "link variable" U^(n) = 
e %a »( n \ related to gauge field by a ll (n) = aeA^fo). 

Monopoles are detected in the lattice using the method 
given by DeGrand and Toussaint [13 : due to the Gauss's 
law, the total magnetic flux emanating from a closed 
surface allows to determine if the surface encloses a 
monopole. The monopole density is defined by half of the 
total number of monopoles and antimonopoles divided by 
the number of elementary cubes in the lattice. We apply 
this definition for both the compact and the noncompact 
formulations of the theory, although some caution should 
be used in this respect. Indeed, monopoles are classical 
solutions of the theory with finite action only for compact 
QED3, where they are known to play a relevant role. In 
the noncompact formulation of QED3 they are not clas- 
sical solutions, but they could give a contribution to the 
Feynman path integral owing to the periodic structure of 
the fermionic sector 1381. 



III. COMPACT VERSUS NONCOMPACT 
FORMULATION 

In order to investigate the onset of the continuum 
physics, it is convenient to consider a dimensionless ob- 
servable and to evaluate it from the lattice for increasing 
(3 until it reaches a plateau. Such an observable can be 
taken to be (3 2 (xx)i which is expected to become con- 
stant in the continuum (/3 —> 00) limit |27l l39| . Numer- 
ical simulations show two regimes: for f3 larger than a 
certain value, the theory is in the continuum limit (flat 
dependence of a dimensionless observable from 0) , other- 
wise the system is in a phase with finite lattice spacing. 
In the former regime, the theory describes continuum 
physics, in the latter one it is appropriate to describe a 
lattice condensed-matter-like system. 

There are a couple of papers by Fiebig and Woloshyn 
in which the two formulations are compared in the finite 
lattice regime 32, 33]. In these papers the /3-dependence 
of the chiral condensate and of the monopole density for 
lattice QED3 with Nf = and Nf = 2 are analyzed for 
both compact and noncompact formulations in the finite 
lattice regime. 

It is shown there that, when (xx) is plotted versus the 
monopole density p m , data points for both theories fall 
on the same curve to a good a ppr oximation (see Fig. E}. 
This led the authors of Refs. |32l l33| to the conclusion 
that the physics of the chiral symmetry breaking is the 
same in the two theories. 

Our program is to study if the conclusion reached by 
Fiebig and Woloshyn can be extended to the continuum 
limit, by looking at the same observables they considered: 
namely the chiral condensate and the monopole density. 
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FIG. 2: Correlation between (xx) an d pm f° r the compact 
(circles) and the noncompact (boxes) theories for Nf = 2 and 
8 3 lattice according to Ref. 



IV. 



NUMERICAL RESULTS 



Since QED 3 is a super-renormalizable theory, the cou- 
pling constant does not display any lattice space depen- 
dence. The continuum limit is approached by merely 
sending j3 = l/(e 2 a) to infinity. In this limit all physical 
quantities can be expressed in units of the scale set by 
the coupling e. Therefore, it is natural to work in terms 
of dimensionless variables such as (3m, L/(3 or /3 2 (xx)> 
which depend on e (L is the lattice size). 

The signature that the continuum limit is approached 
is that data taken at different (3 should overlap on a single 
curve when plotted in dimensionless units |28j . 

In practice, numerical results will not describe the cor- 
rect physics of the system even in the continuum limit 
because of finite volume effects which are particularly 
significant in our case, due to the presence of a massless 
particle, the photon. In principle one should get rid of 
these effects by taking L/j3 — » oo. In practice, this ratio 
is taken to be large, but finite. In Ref. |4£j the authors 
conclude that in order to find chiral symmetry breaking 
for Nf = 2 at least a ratio L/(3 « 5 x 10 3 is required. In 
our simulations the largest value for the L/j3 ratio has 
been 20. 

Our Monte Carlo simulation code was based on the 
hybrid updating algorithm, with a microcanonical time 
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FIG. 3: As in Fig. |5] according to our results 



step set to dt — 0.02. We simulated one flavour of 
staggered fermions corresponding to two flavours of 4- 
component fermions. Most simulations were performed 
on a 12 3 lattice, for bare quark mass ranging in the in- 
terval am = 0.01 -j- 0.05. We made refreshments of the 
gauge (pseudofermion) fields every 7 (13) steps of the 
molecular dynamics. In order to reduce autocorrelation 
effects, "measurements" were taken every 50 steps. Data 
were analyzed by the jackknife method combined with 
binning. 
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FIG. 4: /3 2 (xx) versus f3m in the compact formulation. 

As a first step, we have reproduced the results by 
Fiebig and Woloshyn which are shown in Fig. [3 We find 
that also in our case data points from the two formula- 
tions nicely overlap (see Fig.[3J). It should be noticed that 
data of Fig. [5] were obtained using a linear fit with two 
masses (am=0.025, 0.05) whilst those of Fig.|3have been 
obtained by a quadratic fit with four masses (oto=0.02, 
0.03, 0.04, 0.05), nevertheless the conclusion is the same 
in both cases. We have verified that if we perform a linear 
fit on the subset of our data with masses am=0.02 and 
0.05 and on the subset with masses am=0.03 and 0.05, 
our results nicely compare with those plotted in Fig. [21 
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Then, in Fig. 0] we plot data for P 2 (xx) obtained 
in the compact formulation versus Pm. We restrict 
our attention to the subset of /3 values for which data 
points fall approximately on the same curve, which in 
the present case means j3 = 1.9, 2.0, 2.1, corresponding to 
L/p = 6.31,6.00,5.71. A linear fit of these data points 
gives x 2 /d.o.f. ~ 8.4 and the extrapolated value for 
Pm -> turns out to be P 2 (xx) = (1-54 ± 0.25) x 10~ 3 . 
Restricting the sample to the data at P = 2.1, the 
X 2 /d.o.f. lowers to ~ 1.3 and the extrapolated value 
becomes P 2 (xx) = (0.94 ± 0.28) x 10~ 3 , thus showing 
that there is a strong instability in the determination 
of the chiral limit. If instead a quadratic fit is used 
for the points obtained with P = 1.9,2.0,2.1, we get 
P 2 (xx) = (0.91 ± 0.45) x 10~ 3 with x 2 /d.o.f. ~ 8.7. 
Owing to the large uncertainty, this determination turns 
out to be compatible with both the previous ones. 
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FIG. 5: f3 2 {xx) versus /3m in the noncompact formulation. 

In Fig. we plot data for P 2 (xx) obtained in the non- 
compact formulation versus /3m. Following the same 
strategy outlined before, we restrict our analysis to the 
data obtained with P = 0.7, 0.75, 0.8, which correspond 
to L/P= 17.14,16,15. 

If we consider a linear fit of these data and extrapolate 
to Pm -» 0, we get P 2 (xx) = (0-45 ± 0.03) x 10~ 3 with 
% 2 /d.o.f. ~ 17. Performing the fit only on the data 
obtained with P = 0.8, for which a linear fit gives the 
best x 2 /d.o.f. value ~ 16, we obtain the extrapolated 
value P 2 (xx) = (0-66 ± 0.07) x 10~ 3 . Therefore, also 
in the noncompact formulation the chiral extrapolation 
resulting from a linear fit is largely unstable. A quadratic 
fit in this case gives instead a negative value for P 2 (xx) ■ 

The comparison of the extrapolated value for /3 2 (xx) 
in the two formulations is difficult owing to the instabili- 
ties of the fits and to the low reliability of the linear fits, 
as suggested by the large values of the x 2 /d.o.f. Tak- 
ing an optimistic point of view, one could say that the 
extrapolated P 2 (xx) f° r P — 2-1 in the compact formula- 
tion is compatible with the extrapolated value obtained 
in the noncompact formulation for /3 = 0.8. 



It is worth mentioning that our results in the non- 
compact formulation are consistent with known results: 
indeed, if we carry out a linear fit of the data for 
P = 0.6,0.7,0.8 and am=0.02, 0.03, 0.04, 0.05 and ex- 
trapolate, we get P 2 (xx) = (1-30 ± 0.07) x 10~ 3 with 
an admittedly large x 2 /d.o.f. ~ 20, but very much in 
agreement with the value P 2 (xx) = (1-40 ±0.16) x 10~ 3 
obtained in Ref. |35j . 

We stress again that our results are plagued by strong 
finite volume effects, therefore our conclusions on the ex- 
trapolated values of P 2 (xx) are significant only in the 
compact versus noncompact comparison we are inter- 
ested in. We do not even try to draw any conclusion 
from our data on the critical number of the flavours. As 
a matter of fact a recent paper [2j| shows that, if effects 
are carefully monitored and large lattices, up to 50 3 , are 
used, it is possible to establish that P 2 (xx) < 5 x 10~ 5 . 
For the comparison between compact and noncompact 
QED 3 it is pertinent to carry out the numerical analysis 
with an (approximately) constant value of the ratio L //3. 
This condition is indeed verified even if we performed 
simulations on lattices with fixed (L = 12) size, since the 
range of allowed values for p corresponding to the contin- 
uum limit is narrow (/3 = 1.8 -j- 2.2 in the compact case, 
P = 0.6 -7- 0.9 in the noncompact case). Finite volume 
effects play a "second order" role in our work, since they 
probably only affect the extension of the continuum limit 
window of /3 values. 
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FIG. 6: (xx) versus p m in both the compact and the non- 
compact formulations on a 12 3 lattice. 

In Fig. we plot (xx) versus the monopole density 
p m . Differently from Figs. 12131 it is not evident with the 
present results that the two formulations are equivalent 
also in the continuum limit, although such an equivalence 
cannot yet be excluded. 

In Fig. [7| we plot again (xx) versus the monopole den- 
sity p m , but now on a 32 3 lattice. In this case the chiral 
condensate is extrapolated to zero mass by a quadratic 
fit. In spite of the negative value taken by (xx) f° r large 
P, in this case data for both formulations seem to fall on 
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compact formulations on a 32 3 lattice. 
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In Fig. [8] and Fig. we plot f3 3 p m versus (3m for the 
two formulations; the former quantity is dimensionless, 
therefore, in analogy with the previous cases, we expect 
that data at different (3 values should fall on a single 
curve in the continuum limit. Our results show that this 
is not the case, this suggesting that the continuum limit 
has not been reached for the monopole density. 

Simulations on the 32 3 lattice give practically the same 
results for (3 3 p m , indicating that this observable, unlike 
/3 2 (xx), is volume independent. 

It is important to observe, however, that the monopole 
density is independent of the fermion mass. Since the 
mechanism of confinement in the theory with infinitely 
massive fermions, i.e. in the pure gauge theory, is 
based on monopoles and since the monopole density is 
not affected by the fermion mass, we may conjecture 



that this same mechanism holds also in the chiral 
limit. This supports the arguments by Herbut about 
the confinement in the presence of massless fermion 0,0. 
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FIG. 9: P 3 p m versus f3m in the noncompact formulation. 



V. CONCLUSIONS 

In this paper we have compared the compact and the 
noncompact formulations of QED3 by looking at the be- 
havior of the chiral condensate and the monopole density. 

Numerical results for /3 2 (xx) are compatible with those 
obtained by other groups, although it is still questionable 
if the continuum limit has been reached and if the chiral 
limit is stable. The biggest difficulty for this observable 
is that the chiral extrapolation is rough when a linear fit 
is performed, but gives a negative value when instead a 
quadratic fit is considered. Massive calculations on larger 
lattices are needed to further reduce the finite volume 
effects and to stabilize the chiral limit. 

As far as monopoles are concerned, they appear in 
smaller and smaller numbers for large j3, this making the 
determination of the continuum limit for P 3 p m rather 
problematic. Our results show, however, a very weak 
volume dependence. 

We have analyzed also the relationship between the 
monopole density and the fermion mass, both in compact 
and noncompact QED3. The weak dependence observed 
leads us to conclude that the Polyakov mechanism for 
confinement holds not only in the pure gauge theory, but 
also in presence of massless fermions. 

Finally, we note that, although the chiral condensate 
and monopole density approach the continuum limit in 
two different ranges of /?, the analysis a la "Fiebig and 
Woloshin" does not allow to exclude the equivalence 
of the compact and noncompact lattice formulations of 
QED 3 . 
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